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£N| . Abstract 

q ■ The internal symmetry group U(3,l) of the neutral vector fields with two spins 

[ and 1 is investigated. Massless fields correspond to the generalized Maxwell equa- 

| tions with the gradient term. The symmetry transformations in the coordinate 

space are integro-differential transformations. Using the method of the Hamilto- 
nian formalism the conservation tensors are found, and the quantized theory is 
| studied. The necessity to introduce an indefinite metric is shown. The internal 

symmetry group U(3, 1) being considered, after the transition to electrodynamics, 
reduces to the U(2) group. It is shown that the group of dual transformations is the 
subgroup of the group under consideration. All the linearly independent solutions 
of the equation for a free particle obtained in terms of the projection matrix-dyads. 



o 
o 



9* 1 Canonical formalism 



In the general case, without any constraints, the vector field B^(x) realizes the (0,0) © 
(1/2, 1/2) representation of the Lorentz group and describes four degrees of freedom which 
correspond to states with spins s = and s = 1 (with three spin projections s z = 0, ±1). 
The massive field functions B^(x) satisfy the Klein-Gordon- Fock equation 



3 



X 



(dl - m 2 ) B u (x) = 0, 



where <9 M = (d/dx m ,d/idt, B v (x) = (B n (x) , iB (x)) . 

The corresponding Lagrangian for the neutral fields can be rewritten as follows (within 
unimportant divergent-type terms): 



<~5 



{d,B u f + m 2 Bl . (2) 



The Lagrangian (2) can be connected also with the Stueckelberg formulation of the 
vector field [1]. A Lagrangian of the form (2) also was used [2] in a gauge-invariant 
formulation for a massive neutral vector field. 

Eq.(l) can be represented in the form of first-order equations [3] 

dyip^u - d^ipo + m 2 ^ = 0, 
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*Po = d^, (3) 

with ip^ = B^, ijj = d^B^. In the case m = equations (3) are the generalized Maxwell 
equations with the gradient term (see [4-6]). Now we investigate the symmetry group of 
the four-component neutral vector field B^ which describes two spins 0, 1 (without the 
Lorentz condition, i.e. d^B^ ^ (see Eq.(l)) with the Lagrangian (2) [7]. We call this 
field "the Stueckelberg field" [8]. 

Let us consider the transformations A = (A^) of the field functions B' (x) = A IMU B L ,(x), 
which belong to the group SO(3, 1) but which leave the coordinates x^ unchanged, i.e. 
x 1 = x^. This is different case from that of the Lorentz group where the coordinates x^ 
are transformed. It is easy to verify that the Lagrangian (2) is invariant under this group 
of symmetry transformations as A MQ A M/ g = 5 a p and B'^(x) = B^(x). In accordance with 
the Noether theorem [9,10] the invariance of the action integral under the group of the 
transformations under consideration one yields the law of conservation of the following 
antisymmetric tensor 

Sfj,[aP] = Bpd^Ba - Bad^Bp. (4) 

This tensor coincides with the density of spin momentum corresponding to the con- 
servation law of spin momentum (4): d^S^a^ = 0. 

The invariance of the action integral (2) under the Lorentz transformations of the 
coordinates x'^ = L^ u x u induces transformations of the field functions B' (x 1 ) = L flu B u (x) 
which lead to a conservation law for angular momentum M^^j [9,10], this being the sum 
of the orbital, T^ a x^ and spin, S^a/j], momenta: 

Tfxv = — (d^B a ) {d v B a ) — CSnv, 

= T^aXfl + S^]. (5) 

Therefore we also have here the law of conservation of orbital angular momentum: 
dpT^aXp] = 0. It should be noted that the Lorentz calibration d^B^ = is not invariance 
under our group of symmetry (A^). As the Lorentz calibration extracts the pure spin 1 
of particles we come to the conclusion that the laws of conservation of the orbital, T^xm, 
and spin, angular momenta separately: d^T^xp] = 0, d^S^ap] = are due to the 

presence of two spins and 1; i.e. multi-spin 0, 1. 

Now we consider a wider group of symmetry using the method of the Hamiltonian 
formalism. 

The generalized coordinates in this scheme are 

q,(x) = B,{x). (6) 
The density of the momenta found from Eq.(2) is given by 

M*) = -^TT =% (*)> (7) 
d <? M (x) 
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where 1^ (x) = dq il {x)/dt. The density of the Hamiltonian is defined by the relation- 
ship 

H(x)=ir lt (x)q li (x)-C(x). (8) 

Inserting the density of the Lagrangian (2) and the densities of the momenta (7) into 
Eq.(8) we arrive at 



H{x) = - 



. 2 



(d m B^x) f+Bl (x)+m 2 Bl(x) 
In momentum space the real fields B^(x) are given by 

B,(x) = L- 3 / 2 ]T [B,(k, t)e lk * + fl+(k, t)e 



ikx 



(9) 



(10) 



where L is the normalizing length so that the energy of a quantum is ko = 2n/L and 
the normalizing volume is V = L 3 ; B+(k,t) is the Hermitian conjugated quantity. The 
time dependence of fields in the momentum space is 



fl„(k,t) ~e 



—ikot 



B+(k,t)~ e 



ikot 



(ii) 



where k$ = k 2 +m 2 . Taking into account Eqs.(9), (10) the Hamiltonian H = J H(x)d 3 x 
takes the form 

= £2fcgB /1 (k,t)B+(M). (12) 

k 

Introducing the canonical variables in the momentum space 

g M (k,f) = £ M (M) + £j(M), 

7r M (k, i) =q, (k, i) = -ifo [B„(k, t) - B+(k, t)} , (13) 
the Hamiltonian (12) is rewritten as 



~E [Mk, t)f + kl (9„(k, t)) 2 ] = £ 2fc W, 



(14) 



V? = 



( #i(M) \ 
£ 2 (M) 

£ 3 (M) 
V i£o(M) y 



V?+ = ( Bt(k,t) B}(k,t) Bt(k,t) iB+(k,t)). (15) 



It is obvious that the Hamiltonian (14) is invariant (see [11, 12]]) under the transfor- 
mations of the pseudounitary group U(3, 1) (for real fields B m (k,t), B (k,t)): 



up' = U(p, 



(16) 
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where the complex 4 x 4-matrix U obeys the equation U + U = 1. The matrix of the 
infinitesimal transformations (16) can be represented as 



by 



U — (1 — iuj ) J 4 , +w\p, v ]I\ f j, 1 ,] — ^(jj, v )I(jtv), (17) 
where the antisymmetric i^j and symmetric I^ v ) generators of the group are given 



1/,/* 



2 



(18) 



Here J 4 is the unit 4 x 4-matrix and the elements of the entire algebra e^ u satisfy the 
relations [9] 



,H,v\ 



(19) 



'a/3 ~ ^iwflv^i 

where indexes [/,,v,a,/3 = 1,2,3,4. The parameters of the U(3.1) group obey the equa- 
tions Uq = ^0, ^[ab] — U[abh W *a4] = ~ UJ [aA] , ^* a b) = W H)i W (a4) = ~^(aA) (a,b = 1,2,3). 

The generators of the Z7 (3.1) group: 1I4, I\p,v\, I(pu) have the correct commutation rela- 
tions [13] as may be verified using Eqs.(19). The I\p,v] are the generators of the subal- 
gebra corresponding to the 5*0(3,1) group and il± corresponds to the subgroup U{\). 
Transformations (16) act in the space of the functions (15) and the momentum k is not 
transformed. 

To get the integral of motion (corresponding to the group of symmetry (16) in the 
framework of this formalism one needs to consider the canonical transformations which 
leave the Hamiltonian (14) invariant. Taking into account Eqs.(17), (18), infinitesimal 
transformations (16) can be cast into the form 



1 



-2i 



-UJ{aa)&iJ,v 



B v (k,t), 



SB+(k,t) = iu B+(k,t) + 2u [H B+(k,t)+ 
1 



+ 2i 



B+(k,i). 



(20) 



(21) 



With the help of Eq.(13), transformations (20), (21) can be rewritten in the canonical 
form 

%(k,t) = — vr M (k,t) + 2u [tiu] q l/ (k,t) + 



n u (k,t), 



5vr M (k,t) = -kouoq^t) + 2a; [H 7T l/ (k, *)- 
1 

t 



2k 



g„(k,f). 



(22) 



(23) 
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To get integrals of motion we use the method of generating functions [14]. It is verified 
that the generating function corresponding to transformations (22), (23) is 



F = E{fc(MK(M) + y 



"" (k ' /// ' +M<7m(M)) 2 



ko 



+ 



+^V] [vr^(k,t)g I ,(k,t) - <(k, t)q^(k, i) 



+ 



7^(k,tK(k,t) 
fc 



+ k q^(k,t)q„(k,t) 



(24) 



so that g£(M) = dF/dn'^(k,t), n^(k,t) = dF/dq^(k,t). From Eq.(24) we find the 
conservation law of the following tensors 

J [H = E Mk,t)q u (k,t) - 7i u (k,t)q^(k,t)] = 



Ji 



(/«/) 



= <E[WMk)-C(k)Mt) 



El iu + k oQAK t)q u {k, t)- 



As 



k Q 
A; 



+ fc (g a (k,t)) z 



= E k(k)fe,(k) + KMbM - -5, u (&+(k)6 Q (k 



J 



E 



+ k (q a (k,t)y 



(25) 



(26) 



(27) 



where variables 6 M (k) = y / 2A^_B M (k) and 6+(k) = v / 2A^5+(k) in the quantized theory 
are annihilation and creation operators, respectively. The conserved variables (25)-(27) 
satisfy the equations { J [/H , if} = { J^, if} = {J, H} = 0, where H is given by (14) and 
the {...} are the classical Poisson brackets. In the quantum case, operators 6^(k), b+(k) 
obey the commutation relation [6 M (k), b+(k')] = 5^ u 5 (k — k') and the J. ., J, J are 



the generators of the group of internal symmetry U(3, 1) so that J [fiv] , H 
[J, H] = 0. Using the expansion 

B^x) = L~ 3/2 E [B„(k)e i{k *- kot) + 5+(k) e -* (kx - fcoi) 



(H ' 
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it is easy to verify that the integral of motion, J (Eq.(25)), coincides with the 
conserved spin momentum tensor: 



= 1 J S^d'x. (28) 

So the group SO(3, 1) under consideration (A) with the generators J,, is the subgroup 
of the general group U(3, 1). The conserved quantity J is the number of quanta of the 
field. The generator J defines the subgroup U(l) of the phase transformations and the 
generators J^ v] , J [pv) correspond to the group SU(3, 1). In the general case, transforma- 
tions (16) in the coordinate space are integro-differential transformations and therefore 
the Lagrangian formalism is not convenient for studying this symmetry. The use of the 
canonical formalism and the method of generating functions allow us to investigate in a 
simple manner the group of "addition" symmetry of the vector field which possesses two 
spin values 0, 1 (massive and massless fields). The analogous group of symmetry SU(n) 
of the n-dimensional oscillator has been investigated [11, 12]. The analysis performed is 
readily generalized to arbitrary vector fields. 

So the integral of motion (25)- (27) corresponds to the transformations of internal- 
symmetry (16) which are not induced by the space-time transformations. 

The canonical variables q^(x), n^x) satisfy the commutation relation 

{q^x), n v (x')} t=t , = 6^6 (x - x') . (29) 
Using Eq.(7) we arrive at the relationships 

[q»(x), q v = 6^,6 (x - x') . (30) 

To quantize the fields we should transfer to quantum commutators in accordance with 
the procedure {., -} t=t , — > •], where [q,ir] = qir — nq. Taking into account Eqs.(6), (30) 
we find the commutation relations 

B^x) : B v (x')} =^5(x-x'), (31) 
lt=t' 

It is obvious that commutators (31) correspond to Bose-Einstein statistics. 



2 Quantized fields and indefinite metric 

Now we will consider the quantized theory of the fields with multi-spin 0,1. In accor- 
dance with the general rule that the translation generator P M and the generator of four- 
dimensional rotations M pa are given by [9, 10] 

P il = i J T Ap d 3 x, 
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-i J M 4[pa] d 3 x. (32) 



Using Eqs.(5) and commutation relations (31), it is not difficult to check that genera- 
tors (32) satisfy the commutation relations of the Poincare group. 

Let us discuss the eigenvalues of the operator of the energy P = —iP 4 of the fields 
considered. In the momentum representation, the energy operator has the form: 

P = / h (&+(k)6„(k) - & + (k)6 (k)) d% (33) 

where operators 

6j(k) = ^2k B+(k), 6„(k) = v^^(k) (£„(k,t) = £? M (k) exp(-iM)) , 
as follows from Eqs.(31), must obey the following commutation relations: 

6 M (k),6+(k')] = V<J(k-k'). (34) 

Commutation relation (34) is invariant under the group of symmetry (16) because the 
operator of finite transformations U is unitary. It is seen from Eq.(33) that the field 
energy, in classical theory, is not positive-definite. In commutation relations (34) there is 
a minus sign (as 64 (k) = ib (k), 64 (k) = ib^ (k)), and the operators bo (k), b£ (k) obey the 
"incorrect" commutation relation 

& (k),&+(k')] =-5(k-k'). 

From this equation it is seen that the term which appears in the energy operator with 
a (— ) sign satisfies the "incorrect" commutation relation. The "incorrect" commutation 
relations are not compatible, however, with the assumption that the fields are real. As 
in the case of the electromagnetic field, this difficulty is surmounted by introducing an 
indefinite metric (see e.g., [15]). There are two possibilities. According to the first one we 
can consider operators &o(k), b+ (k'), as the creation operators, and b a (k), b^ (k'), as the 
annihilation operators of particles. The vacuum state | 0)i is defined by the requirement: 

6 (k) I 0)i = 6+(k) I 0)i = 0, 

i(0 I 0)i = I. (35) 
The basis in Hilbert space is given by 

I m,n>i = -7^ ():(k))>o(k)r I OH (36) 

with the normalization condition 

(m, n I m', n') 1 = 5 mm ,5 nn >. (37) 
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In this case there is Hilbert space but the eigenstates of the energy operator (33) are not 
positive defined and the interpretation of the states with negative energy is problematical. 

The second possibility [15] is more favorable and connected with introducing an indef- 
inite metric. In this case the algebra (34) can be represented as the operator algebra in 
the states with the indefinite metric and here operators b+(k) are the creation and &^(k) 
the annihilation operators. As the basis of the irreducible representation in Hilbert space 
we choose the following vectors: 



m, n) = 



^(«(k))">+(k))"|0>, 
ym\n\ 



6„(k) | 0) = 0, (0 | 0) = 1. (38) 

With the help of Eqs.(34) we arrive at 

(m,n \ m',n') = (-l) n 8 mm/ 5 nn >. (39) 



In this way, the state space of vector fields has an indefinite metric, i.e. the square 
vector norm can be negative. Using Eqs.(34), it is easy to verify that the eigenvalues of 
the operator P are positive; however, the metric is still indefinite. 

As usual, in the theory with an indefinite metric, it is necessary to divide the space 
into "physical" and "nonphysical" subspaces (see [15, 16]). The "physical" subspace cor- 
responds to a positive square norm, and the "nonphysical" to a negative square norm. 
The state vectors | to, 0) create the "physical" subspace H p and | to, n) (at n^0) create 
the "nonphysical" subspace H n . The physical states permit the usual probability inter- 
pretation. It is apparent from Eq.(39) that the "physical" and "nonphysical" subspaces 
are orthogonal (to, | to', n') = 0, and therefore the state space is the direct sum of the 
two subspaces H p and H n . We can then represent any vector in the form 

I) HU |>n, (40) 

where |) p G H p , \) n G H n . In scattering processes, states | in), \ out) should belong 
to the "physical" subspace H p . It is possible also to consider transitions between states 
which belong to the "nonphysical" subspace H n . But the consideration of transitions 
between the vectors \) p , |)„ presents difficulties [15] because such transitions violate the 
unitarity of the ^-matrix and the usual probability interpretation. 



3 Internal symmetry of the electromagnetic field 

It will be shown that the symmetry-group of an electromagnetic field is U(2). When 
we impose the constraints d^B^ = (The Lorentz condition), to = 0, we arrive at the 
transition to the formulation of electrodynamics, and the Hamiltonian (33) takes the form 

H = Y, k o (tf(k)&i(k) + bt(k)b 2 (k)) . (41) 

k 
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In this case the Hamiltonian (41) is invariant under the transformations of the U (2) 
group: 

raHfi— DCS!)' 

where n 2 = 1, a, 9 are real group parameters, r are the Pauli matrices. So, the internal 
symmetry group U(3, 1) being considered, after the transition to electrodynamics, reduces 
to the U(2) group. Using the procedure described above, after finding the generalizing 
function, we come to the following conserved quantities: 

J 1 = V[6+(k)6 2 (k) + 6+(k)6 1 (k)" , 
z k 

J 2 ^£[^(k)Mk)-^(k)6 2 (k)" , 
z k 

J 3 = ^K( k )Mk)-^(k)6 2 (k)], (43) 
z k 

J o = lH K(k)&i(k) + 6 2 + (k)6 2 (k)] . (44) 
z k 

The generators J, satisfy the commutation relation of the SU (2) group, and com- 
mute with J : [Jj, Jo] = 0. There is a representation of the rotation generators in the 
Hilbert space in the form of Eqs.(43) in [17, 18]. The average values of the operators 
(43) are identified in [19] with the Stokes parameters characterizing the polarization of 
the electromagnetic wave. 

It is obvious, that the transformations (42), and the integrals of motion (43), (44), can 
be written in coordinate space but in nonlocal form. So, in the coordinate representation 
the expression (44) takes the nonlocal form [20] 

Iff E(x)E(y) + H(x)H(y) 

The transformations (42) at 9 = 0, which lead to the integral of motion (45) are local 
only in the momentum representation, and in the coordinate representation are integro- 
differential transformations. 

It should be noted that the transformations (42) with rii = n 3 = 0, a = 0, n 2 = 1, 
lead to the dual transformations of the strengths of electric E, and magnetic H fields: 

6 9 
E' = Ecos- + Hsin-, 
2 2 

H' = H cos - — E sin - . (46) 
2 2 v ' 

Hence, the well known group of dual transformations is the subgroup of the SU(2) 

group under consideration. 
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4 First-order equations and density matrix 

Now we consider the matrix formulation of the first-order of the Stueckelberg fields (3) 
which is convenient for constructing the density matrix and for some electrodynamics 
calculations. All the linearly independent solutions of the equation for a free particle will 
be obtained in terms of the projection matrix-dyads. 
Let us introduce the 11-dimensional function 



I ( -^o 



m 



m% j (A = 0,aUH), (47) 



where v — 1, 2, 3, 4. Using the elements of the entire algebra Eq. (19), equations 
(3) can be written in the form of one equation 



£ i*,n _|_ £ o,o _|_ 1 £ [H.[H 



+ m 

After introducing 11-dimensional matrices 

1 



*b(x) = 0. (48) 

AB 



1 1 1 ;=fc -^ + e o,o + _ e [H,[H ) (49) 
2 



Eq. (48) takes the form of the relativistic wave equation of the first order: 

(a^ + m) = 0. (50) 

We took into account that in in Eq. (49) is the unit matrix in 11— dimensional space. 

Eq. (50) represents the Stueckelberg equation for massive fields in the matrix form. 
When fields ^ a(x) are complex values, Eq. (50) describes charged particles with multi- 
spin 0, 1. 

It should be noted that the matrices a M can be represented as 

^ = ^ ] +^\ 

0(1) _ £ m,[H _|_ £ [H>m 

/3(o) = e ".o + (51) 

where the 10— dimensional @jp and 5— dimensional (3^ matrices obey the Petiau- 
Duffin-Kemmer [21-23] algebra: 

P^PvPa + PaPvPv = S^(3 a + 5 au (3^ (52) 
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so that the equations for spin-1 and spin-0 particles are 



+ m) = 0, = 1 ( ^ ) , (53) 



+ m)*W(x)=0, * ( °^) = ^( m ^V (54) 



The 10— dimensional Petiau-Duffin-Kemmer equation (53) is equivalent to the Proca 
equations [24] for spin-1 particles and the 5— dimensional Eq. (54) is equivalent to the 
Klein-Gordon- Fock equation for scalar particles. The 11— dimensional Eq. (50) de- 
scribes fields with two spins 0, 1. It is not difficult to verify (using Eqs. (19)) that 
the 11— dimensional matrices (49) satisfy the algebra (see also [25]): 

= 2 (Sfj, v a a + 5 av a^ + ^a,) . (55) 

This algebra is more complicated than the Petiau-Duffin-Kemmer algebra (52). Differ- 
ent representations of the Petiau-Duffin-Kemmer algebra (52) were considered in [26-28]. 

Now we find the solutions to Eq. (50) corresponding to definite values of the energy 
and momentum of a quantum of the massive fields. In the momentum space Eq. (50) 
becomes 

— ip^p = em^ p , (56) 

where p = a^p^, p^ = (p,ipo), p 2 = p 2 — pi = — m 2 ; the value of e — 1 corresponds 
to positive energy and e — — 1 to negative energy. Here p means the momentum of a 
field-quantum. It may be verified using (55) that the equality 

p? = p 2 p (57) 

is valid. Following the general method of projection operators [29, 30], we find solutions 
to Eq. (56) in the form of the projection matrix 

M, - ^f^, (58) 

so that 

M 2 e = M £ , (59) 

and e — ±1. Every column of the matrix M £ can be considered as an eigenvector ^ p of 
equation (56) with eigenvalue em. Eq. (59) for projection operators tells that matrix M £ 
can be transformed into diagonal form, with the diagonal containing only ones and zeroes. 
So the M £ acting on the wave function \1/ will retain components which correspond to the 
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eigenvalue em. The generators of the Lorentz group in the 11— dimensional space being 
considered are given by 

J^ = )9fW ) -)9W ) . (60) 

It should be noted that matrices (60) act in the 10— dimensional subspace (mip^, Vvh) 
because the scalar if) Q is an invariant of the Lorentz transformations. So matrices (60) 
are also generators of the Lorentz group for the Petiau-Duffin-Kemmer fields of Eq. (53). 
Using properties (19), we get the commutation relations 

\Jpai J pv\ ^ap'Jpu ^pu'Jap $ ppJ ay ^ai/'Jppi (^1) 

[ax, Jpu] = hp^u - huCtp. (62) 

Relationship (61) is a well known commutation relation for generators of the Lorentz 
group 5*0(3, 1). Equation (50) is form-invariant under the Lorentz transformations since 
relation (62) is valid. To guarantee the existence of a relativistically invariant bilinear 
form 

^ = y+rjty, (63) 

where is the Hermitian-conjugate wave function, we should construct a Hermitian- 
izing matrix rj with the properties [9, 28, 30]: 

rjai — —aiT), r\a^ = a.tf) (i = 1,2,3). (64) 

Such a matrix exists and is given by 

„ _ c 0,0 , 9/ q(1)2 t 

rj — —e + zp 4 — iio, 

/ 10 = £ ^ + i £ [H,[H ! ( 65 ) 

where the matrix rf^ = 2(3^ 2 — Iiq plays the role of a Hermitianizing matrix for 
the Petiau-Duffin-Kemmer equation (53) [9]. The operator of the squared spin (squared 
Pauli-Lubanski vector) is given by 

/ 1 \ 2 1 

° 2 = \ ^; £ »™PPv J <*p) = 7^2 ( J lvP 2 ~~ J P^ J vaPpPv) • (66) 



It may be verified that this operator obeys the minimal equation 

a 2 (a 2 - 2) = 0, (67) 

so that eigenvalues of the squared spin operator a 2 are s(s + 1) = and s(s + 1) = 2. 
This confirms that the considered fields describe the superposition of two spins s = and 
s — 1. To separate these states we use the projection operators 

S?o) = 1 - y , S 2 (l) = £ (68) 
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having the properties S^S 2 ^ 



0, (5?- < 



(o)J 



C2 



C2 
3(1) 



C2 o2 _i_ Q2 
°(1)> D (0) + D (1 



(1) 



1, where 



1 = In is the unit matrix in 11— dimensional space. In accordance with the general 
properties of the projection operators, the matrices S? \, S?^ acting on the wave function 
extract pure states with spin and 1, respectively. Now we introduce the operator of the 
spin projection on the direction of the momentum p : 



2 | P 



f ^abcPa J be 



TtabcP. 



•AM*, 



(69) 



where | 
operator is 



Pi + Pi + Pi- The minimal matrix equation for the spin projection 



Op(op - 1) (a p + 1) = 
and the corresponding projection operators are given by 
* 1 



'(±i) 



-a p [a p ±l) 



'(0) 



I -at 



Operators £(±1) correspond to the spin projections s p = ±1 and 5(0) to s f 



(70) 

(71) 
0. It is 



easy to verify that the required commutation relations hold: 



S, 



(0)'P 



o. 



(i),P 



r-2 



S(o),P 
(0)> ^(o) 



= 0, 
0. 



(72) 



^(0))^(±1)J - P(l),^(±l) 

Thus the projection matrices extracting pure states with definite spin, spin projection 
and energy take the form 



A (i) 



ip{ip-em)\ 

2m 2 2°> {ap ± 1} ' 



ip (ip — em) a 1 



A (o) 



2m 2 
ip (ip — em) 
2m~ 2 



1 - 



(73) 



where we took into account that (a 2 /2)a p 



a p . 



Projection operators A £±1 , A' 1 ' 



extract states with spin 1 and spin projections ±1, 0, and corresponds to spin 0. 
The A £j ±i, AW, A^°) are the density matrices for pure spin spates. It is easy to consider 
impure states by summation of Eqs. (73) over spin projections and spins. Projection 
operators for pure states can be represented as matrix-dyads [29, 30]: 



A £ ,±i = ^ 



e,±l 



■ ^ 



£,±1, 



Af» = V £ ■ 



(74) 



where the wave functions ^,±1, are the solution of the field equations for spin 1 
and spin projections ±1 and 0, respectively, and corresponds to the solution with 
spin 0. Expressions (73), (74) are convenient for calculating different electrodynamics 
processes involving polarized vector charged particles. 
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